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CHAPTER  1 


INTRODUCTION  AND  SUMMARY 

This  final  report  on  AFOSR  Grant  No.  77-3462^  is  a 
technical  description  of  exterior-interior  aperture  coupling 
of  a  plane  electromagnetic  wave  to  a  rectangular  cavity  con¬ 
taining  a  wire  obstacle.  Portions  of  the  research  have  been 
described  in  previous  annual  and  interim  reports.  In  addition, 
portions  have  been  presented  at  meetings  of  the  International 
Scientific  Radio  Union  (URSI)  and  in  a  paper  published  in 
Radio  Science  (Johnson,  Howard,  and  Dudley,  1979).  All  of  the 
important  analysis  and  results  for  the  entire  program  are  con¬ 
tained  in  this  final  report.  In  addition,  we  contemplate  one 
additional  journal  article  based  on  the  final  results  contained 
herein. 

The  problem  of  coupling  electromagnetic  energy  into  a 
cavity  through  an  aperture  and  excitation  of  current  on  a  wire 
in  the  interior  is  one  requiring  considerable  care  both  in 
formulation  and  numerical  reduction.  Essentially  the  procedure 
we  have  used  is  to  begin  with  Maxwell's  equations,  produce 
vector  Helmholtz  equations  for  the  electric  and  magnetic  fields, 
and  invert  the  ensuing  second  order  vector  differential  operators 
through  application  of  a  dyadic  Green's  theorem.  The  results 
are  expressions  for  the  electric  and  magnetic  fields  that  can 
be  specialized  as  integral  equations  by  taking  the  observation 
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points  to  the  aperture  and  to  the  surface  of  the  interior  wire. 
The  coupled  integral  equations  produced  are  next  solved  by  the 
method  of  moments,  a  procedure  which  transforms  the  integral 
equations  into  a  matrix  equation  that  can  be  numerically 
inverted . 

To  simplify  the  technical  bookkeeping,  we  have  formu¬ 
lated  the  problem  in  dyadic  notation.  This  approach  is 
particularly  useful  in  displaying  the  various  components  in  the 
Green's  functions.  One  of  the  significant  contributions  of  this 
study  has  been  in  the  careful  construction  of  the  dyads  and 
their  evaluation  in  regions  where  they  become  singular.  This 
subject  is  treated  in  detail  in  our  journal  article  recently 
published  in  Radio  Science,  referred  to  above.  We  have  been 
particularly  careful  concerning  the  distributional  character 
of  some  of  the  singularities  and  the  completeness  of  the  Green's 
dyadic  eigenfunction  expansions,  a  fact  that  has  led  to  some 
controversy  in  the  literature.  We  feel  that  our  contributions 
have  resolved  some  confusing  statements  on  the  longitudinal 
wave  functions  and  we  are  confident  of  our  results. 

The  numerical  reduction  has  occupied  a  major  portion  of 
our  time  during  the  last  year  of  the  research.  We  have  develop¬ 
ed  a  summation  method  for  evaluation  of  otherwise  slowly 
converging  eignefunction  series  and  thereby  have  been  able  to 
significantly  reduce  the  numerical  processing  time.  We  have 
discovered  a  sign  discrepancy  between  our  work  and  that  of  Seidel 
(1978) .  Our  analysis  supports  the  sign  obtained  in  our  results. 


In  Chapter  2,  we  discuss  the  problem  geometry  and 
produce  the  formulation  in  the  frequency  domain.  In  Chapter  3, 
we  discuss  the  scheme  for  numerical  solution  with  particular 
attention  to  subdividing  the  matrix  to  be  inverted  by  the 
method  of  moments.  In  Chapter  4,  we  discuss  the  special  tech¬ 
niques  developed  to  sum  the  series  appearing  in  the  various 
elements  of  the  matrix.  In  Chapter  5,  we  give  sample  numerical 
results  for  aperture  fields  and  wire  currents.  We  have  included 
an  appendix  describing  our  analysis  of  the  sign  difference 
between  this  work  and  that  of  Seidel. 

As  a  footnote,  we  should  like  to  point  out  a  dissapoint- 
ment  in  one  phase  of  the  work.  We  had  hoped  to  be  able  to 
provide  some  transient  results  for  the  wire  current  by  obtaining 
the  currents  at  various  positions  along  the  wire  over  a  wide 
range  of  frequencies,  multiplying  by  the  spectrum  of  a  suitable 
input  pulse,  and  taking  the  numerical  inverse  Fourier  transform. 
Indeed,  Figures  21  and  22  are  representative  of  our  efforts  in 
this  regard.  Unfortunately  the  resonances  are  so  narrow  (in 
engineering  terms,  the  Q  is  extremely  high)  that  we  could  not 
obtain  a  causal  inverse  transform.  In  other  words,  the  resonance 
is  so  localized,  it  is  as  if  we  were  attempting  to  take  the 
inverse  transform  of  a  delta  function.  We  verified  our  difficul¬ 
ties  by  creating  an  RLC  circuit  model,  solving  for  the  current 
in  the  frequency  domain,  and  attempting  a  numerical  inverse  trans 
form.  All  goes  well  until  one  selects  a  value  for  the  resistance 
small  enough  to  again  produce  a  very  high  Q  situation.  At  this 
point  all  codes  we  have  tried  for  performance  of  the  numerical 
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inverse  start  failing.  We  believe  that  this  problem  can  be 
solved  by  analytic  continuation  of  the  frequency  variable  oj 
away  from  the  real  axis  to  avoid  the  poles  of  the  transfer 
function  that  are  located  close  to  the  axis.  Since  attempting 
such  a  procedure  requires  a  major  set  of  changes  in  our  computer 
program  we  have  not  tried  to  do  this. 


CHAPTER  2 


FREQUENCY  DOMAIN  FORMULATION 

The  frequency  domain  version  of  the  boundary  value 
problem  of  Figure  1  is  now  considered  by  specializing  the 
exterior  source  to  that  of  a  time  harmonic  elut  plane  wave.  The 
problem  formulation  of  this  section  will  result  in  a  numerically 
tractable  system  of  coupled  integro-dif ferential  equations  for 
the  thin  wire  currents  and  tangential  aperture  electric  fields. 

Once  these  wire  currents  and  tangential  aperture  electric  fields 
are  determined,  fields  anywhere  in  the  exterior  or  interior  regions 
may  be  computed  by  use  of  the  appropriate  Green's  functions. 

First  potential  dyads  for  cavity  interior  and  exterior 
half-space  are  constructed.  Next  these  dyads  are  used  to  obtain 
integral  representations  of  the  electromagnetic  fields  in  each 
region.  These  integral  representations  are  given  in  terms  of 
the  unknown  tangential  aperture  electric  fields  and  wire  currents. 
Continuity  of  the  aperture's  tangential  electric  fields  and  the 
boundary  condition  that  the  tangential  component  of  the  electric 
field  is  zero  along  the  wire  are  enforced  to  yield  a  system  of 
coupled  integro-dif ferential  equations  for  the  unknown  aperture 
electric  fields  and  wire  currents. 

The  potential  Green's  dyads  for  the  exterior  half-space 
and  cavity  interior  are  now  constructed.  The  frequency  domain 
solution  of  Maxwell's  equations  driven  by  electric  currents  may 
by  expressed  as 
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where 


E  = 


1 

iwey 


2  2 
-  (Vz  +  1c 

n  x  A  = 

V  -  A  = 


(k2  +  VV* )  A 


)A  ■  yj 

0  on  a  perfect  electric 
n  conductor  (p.e.c.). 


(2.1a) 

(2.1b) 


(2.1c) 


In  (2.c)  n  denotes  a  unit  normal  vector. 

The  radiation  condition  at  infinity  and  Meixner's  edge  condition 
[Meixner  1972]  are  imposed  where  appropriate  to  insure  a  unique 
solution.  Similarly  the  solution  of  Maxwell's  equations  driven 
by  magnetic  currents  is  given  by 


where 


H  = 


E  = 


lweu 


<*2  ♦ 


VV*)F 


VxF 

e 


-  (V2  +  k2)F  =  eM 

n  x  (VxF)  =  Oi  on  a  p.e.c. 
n  •  F  =  0  1 


(2.2a) 

(2.2b) 


(2.2c) 


Solutions  of  the  above  systems  (2.1)  and  {2.2)  may  be 
obtained  by  use  of  the  corresponding  magnetic  and  electric  vector 
potential  dyads  GA  and  gp  .  The  magnetic  potential  dyad  satisfies 

-(V2  +  k2)  Ga(?|P)  =  I  5  (r-r ' )  (2.3) 

n  x  GA(r|r')  =  o|  on  a  p.e.c. 

V  •  Ga(?|P)  »  0) 
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where  I  is  the  identy  dyad.  Similarly  the  electric  vector  poten¬ 
tial  dyad  satisfies 

-(V2  +  k2)iF(?|?‘)  =  I  S(F-F')  (2.4) 

n  x  (V  x  gF(F|F' ))=  0|  on  a  p.e.c. 
n  *  gF(F|“’ )  =  0  I 

Since  these  potential  dyads  are  diagonal  the  exterior 
half-space  dyads,  denoted  by  the  superscript  "hs" ,  may  be  con¬ 
structed  by  image  theory.  The  results  are  presented  in  Table  1. 
Since  systems  (2-3)  and  (2-4)  determine  self-ad jont  boundary 
value  problems,  the  interior  cavity  dyads,  denoted  by  superscript 
"c" ,  may  be  given  in  terms  of  their  complete  eigenfunction 
expansions.  These  expansions  are  given  in  Table  2.  The  corre¬ 
sponding  field  dyads  may  be  obtained  by  dif ferenation,  however 
they  are  highly  singular  and  should  be  interpreted  as  distributions 
(Johnson,  Howard,  and  Dudley ,  1979) . 


Table  1.  Potential  Dyads  for  The  Exterior  Half-space 


hs 

GA(r|r' ) 


ii;(r|r,)-it»i(r|r' )  0  0 

0  <l>  (r  |r* )  (r  |r ' )  0 

0  0  ^(r|r' )+^i(r|r') 


-hs  _ 
gF(r|r')  = 


!t»(r|r' )+^i(r|r' )  0  0 

0  (r  |  r' )  (r  |r ' )  0 

__  0  0  i|>  (r  |r' )  -4^  (r  |r' ' ) 


where 


and 


-ikR 


4»(r|r') 


4irR 

-ikR. 


r  *  [(x-x')2  +  (y-y' ) 2  +  (z-z')2l 
R.*  [(x-x')2  +  (y-y')2  +  (z+z'-c)2] 


The  problem  geometry  is  given  in  Figure  1. 


Table  2.  Potential  Dyadic  Green's  Functions  for  A  Rectangular 

* 

Cavity 


c  ,  *  e  e 

G  (r  7* )-  — _ ■a— 

Vr  ’  abc  V  2  .2 

Kj2an  "  * 


(cc)x(ss)y(ss)z 


(ss)  x (cc )  (ss)  z 


(3s)  x  (ss)  y  (cc) 


<3  -V  e.  e  £ 

^  (r  r* )  ■  — — S — 2 

abc  £,ra,n=0„  2  2 

2mn  "  * 


(ss)  (cc)  (cc) 
x  y  z 


(cc)  (ss)  (cc) 
x  y  z 


(cc)x(cc)y(ss)z 


where  (cc)  (ss)  (ss)  ■  cos (k  ji  cos(k  x* )  sin(k  y)sin(k  y*)sin(k  z)sin(k  z*),etc 
x  y  2  x  x  y  y  z  z 

,  mir  ,  rnr  £tt  2  .  2  , 2  ,  2  ,  2 

k  -  -  ,  k  -  —  ,  k  «  —  ,  K0  “k  +k  +k  -  k 

x  a  y  b  z  c  Juan  x  y  z 


1,  if  i  *  0 

i  2,  if  i  4  0 


The  cavity  dimensions  are  a,b,  and  c  in  the  x,  y,  and  z 
directions  respectively. 
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Now  that  the  appropriate  Green's  dyads  have  been  con¬ 
structed,  they  will  be  used  to  obtain  integral  representations 
of  the  electromagnetic  fields  in  each  region.  The  equivalence 
principle  [Harrington  1961,  p.  106]  is  used  to  obtain  an 
equivalent  problem  for  the  exterior  half-space  as  illustrated 
in  Figure  2.  In  the  equivalent  problem  the  aperture  has  been 
shorted  and  the  magnetic  surface  current  Mg  has  been  placed  on 
the  surface  where  the  aperture  was  previously  located.  The 
magnetic  surface  current  is  chosen  to  make  the  electric  field 
jump  to  its  origional  value  just  above  the  aperture  surface. 

The  equivalent  problem  for  the  interior  region  is  obtained 
similarly  and  is  illustrated  in  Figure  3. 

It  is  convenient  to  separate  the  electromagnetic  field 
in  the  exterior  half-space  into  two  terms.  The  short  circuit 
fields,  denoted  by  the  superscript  "sc",  are  those  fields  that 
would  exist  if  the  aperture  were  not  present.  These  fields  may 
easily  be  found  by  geometric  optics.  The  fields  due  to  the 
presence  of  the  aperture  or  equivalently  scattered  by  the  mag¬ 
netic  surface  current  may  be  computed  by  use  of  the  electric 
_  hs 

potential  dyad  gp  of  Table  1  and  (2.2).  In  particular  the 
magnetic  field  in  the  exterior  half-space  is  given  by 

hs  _ 

g p ( r | r ' )  Mg(r')  ds'  (2.5) 

A 

where  A  is  the  aperture  surface.  The  first  term  on  the  right 
of  (2.5)  is  the  short  circuit  magnetic  field,  the  second  is  the 


H  (r) 


_sc 

H(r) 


ioju 


(k2+VV« 
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magnetic  field  due  to  the  presence  of  the  aperture.  The  electric 
field  for  the  exterior  half-space  may  be  found  similarly  by  use 
of  (2.2b)  and  (2.2c). 

The  fields  in  the  cavity  interior  (Figure  3)  are  the  sum 
of  fields  due  to  magnetic  surface  current  plus  the  fields  scattered 
by  the  wire.  These  fields  may  be  obtained  by  use  of  the  cavity 
potential  dyads  (Table  2)  and  (2.1)  through  (2.4).  The  magnetic 
field  in  the  cavity  interior  is  given  by 


H (r) “  -^-(k2+VV‘) 
i<uy 


rr 


gpUlrMHgfc'Jds'  +  V  x 


GA(r|r* )  J(r' )ds' .  (2.6) 


The  electric  field  in  the  cavity  interior  may  be  found  in  a 
similar  manner. 

Note  that  by  construction  the  equivalent  problems  (Figures 
2  and  3)  satisfy  continuity  of  the  electric  field  across  the 
aperture.  Before  enforcing  continuity  of  the  aperture's  tangential 
magnetic  field,  it  is  recalled  that  under  the  thin  wire  assumptions 
of  Section  1  all  current  is  directed  along  the  wire  axis.  Thus 
continuity  of  the  apertures  tangential  magnetic  field  plus  (2.5) 
and  (2.6)  yield 


2iuiy  H 


inc 

y 


(r) 


(k2+ 


A 


hs 

tgF( 


ill 


c'yx,_ 

+  gp  (r | r ’ ) ]Ex(r' )ds' 


(2.7a) 
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and 


inc  -  a2 
-2ioiy  H  (r)  =  (k  +  —  *■) 

x  ax^ 


hs^xx 

tgF(r|r*) 


c,xx 

+  gp(r I r ' ) ]  Ey (r' )ds' 


-a2 

3y3x 


A 

hsz.YZ  C'ZY_  _  a 

[gF(rlr,)+gF(r|r,)]Ex(r,)ds,+ia>u^- 


(2.7b) 

C'1Y- 

GA  ( r ]  r ' )  J  Cr ' )  ds '  . 


A  W 

The  tangential  component  of  the  electric  field  along  the  wire  is 
set  to  zero.  This  gives 


0 


1 

idie 


(k2+ 


■f  c,yy 

Ga(?|?,)J  (?)ds' 
W 


c,xx_ 

gp (r | r ’ ) Ey (r ' ) ds ’ (2. 
A 


Equation  set  (2.7)  is  a  coupled  system  of  integro-dif ferential 
equations  for  the  unknown  aperture  fields  and  wire  currents. 

In  the  next  section  an  efficient  numerical  solution  scheme  for 
this  system  will  be  presented. 


7c) 


CHAPTER  3 


FREQUENCY  DOMAIN  SOLUTION  SCHEME 


A  numerical  solution  scheme  for  the  linear  system  of 
coupled  integro-differential  equations  (2.7)  derived  in  the 
previous  section  is  now  presented.  The  method  of  moments 
[Harrington  1968]  with  pulse  expansions  and  point  matching  is 
used.  As  illustrated  in  Figures  4  and  5,  the  unknown, 
tangential,  aperture  electric  field  components  are  approximated 
by  piecewise  constant  rectangular  patches.  The  x  component  of 
the  electric  field  is  expressed  as 

NE 

Vr')  -  £*  exi  Fxi(r,>  <3-8a> 


where 


P  , (r* )= 
xi 


if  r'  is  in  the  i  rectangle  of  Fig.  4 
otherwise  , 


exi  are  unknown  constant  coefficients, 
of  nonzero  Ex  pulses.  Similarly 


Vr’> 


NE 

-  Ly 

i=l 


'yi 


P  . (r1 ) 
yi 


and  NE  is  the  number  of 
x 


(3.8b) 


where 


(r  ’ ) 


1,  if  r1  is  in  the  ith  rectangle  of  Figure  5 
0,  otherwise  , 


e^  are  unknown  coefficients,  and  NE^  is  the  number  of  nonzero 
Ey  pulses. 
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Note  that  the  Ex  and  Ey  aperture  pulse  expansions 
(Figures  4  and  5)  are  offset.  Centers  of  the  Ex  pulses  corre¬ 
spond  to  corners  of  Ey  pulses  and  vice  versa.  This  will  later 
allow  efficient  evaluation  of  the  x  and  y  derivatives  in  (2.7a) 
and  (2.7b)  by  finite  difference  techniques.  Wilton  and  Glission 
[.1976  ]  applied  this  technique,  offsetting  the  expansion  functions, 
to  the  problem  of  scattering  from  a  rectangular  plate. 

Since  Jy  is  assumed  uniform  about  the  wire  circumference 
it  is  helpful  to  define  the  current  I  by 


I(y')  =  2tt  rw  Jy(y') 


(3.9) 


where  rw  is  the  wire  radius.  Under  the  thin  wire  approximation 
all  current  is  assumed  to  run  along  the  wire's  center  [Seidel  1977, 
1978].  This  current  is  now  approximated  by  (Figure  6) 


where 


I  (y' ) 


NW 

=  £ : 

i=l 


Pwi(y,) 


(3.10) 


pwi<Y'> 


1,  if  y'  belongs  to  the  i  interval 
0,  otherwise 


1^  are  the  unknown  current  coefficients,  and  NW  is  the  number 
of  nonzero  wire  current  pulses. 

Enforcement  of  (3.7a)  at  the  nonzero  Ex  aperture  pulse 
centers,  (3.7b)  at  the  nonzero  Ey  aperture  pulse  centers,  and 
(3.7c)  at  the  wire  surface  points  (xc  +  rw,ym,  zc) ,  where  xc  and 
zc  are  the  respective  x  and  z  corrdinates  of  the  wire  center  and 


y  (m  =  1,  NW)  are  the  nonzero  wire  current  pulse  center  points, 
m 

yields  the  NE  +  NE  +  NW  square  matrix  equation 
x  y 


(3.11) 


The  T  submatrices  are  given  by 


1 1  -  .2  [  [  hs _i_yy_  C'ZY,_ 

=  (k2+-^-j)  JJ  [gF(rmlr,)+  gF(rm|r,)]Pxi(r,)dx,dy' (m=l,NEx; 


i=l,  NEX) 


,  -  .2  hs,xx  c,xx  __ 

Tii  *  dby  j  J  ^rnJ  r*  ^  ^nj  r'  ^  ^Pyi  (r '  J  ^x'dy' 

A  NEX?  i=l,  NEy) 


(3.12a) 


(3.12b) 


-a 

3x3y 


’  f  hs^yy  c,yy  _ 

[gF(rm|r')+gF(rm|r,)]Pxi(r')dx,dy' (m-1, 

NEy?  i*l,NEx) 


(3.12c) 


0  a2  hs,xx  c,xx  _ 

=  (k  +-^)  [gF(rm|r' )+gp(rm|r' )Pv.  (r' )dx'ds'  (m=l. 


F  m1  yi 


NEy;i=l,NEy) 


.  C'W 

GA  (Fm!?,)pwi(y,)dy'd^  (m-l,NEy;i=l,NW) 


(3 . 12d) 


(3 . 12e) 


i 

i 

i 
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-  C  ,  XX 

Tmi  =  Tz  jj^F(?mlF,)  Pyi^Mdx'dy'  (m-l,NW; i=l,NEy) 
A 


(3. 12f  I 


T33  =  -r— i- 
mi  line. 


2  s2  ciyy_ 

r-(kz+  -2-)  G(rJr')P  .  (y')dy'd4>' (m=l,NW;i=l,NW)  (3.12g) 
■ '  av  JJ  m  W1 


The  forcing  function's  components  are 


f~  =  2iwy  H  (r  )  (m=l,NEv) 
m  y  m  x 


f^  =  -2iujy  Hx  <rm)  (m=-l,NEy) 


(3.13a) 


(3.13b) 


Once  the  T  matrix  elements  are  obtained  the  unknown  aperture 
field  and  wire  current  coefficients  may  be  calculated  by 
Gaussian  elimination. 

Evaluation  of  the  T  matrix  elements  (3.12a  -  3.12d)  is 
now  considered.  To  avoid  kernels  that  becomt  highly  singluar 
as  approaches  r ' ,  the  derivatives  are  kept  outside  the  integral 
signs  and  are  evaluated  by  finite  differences.  The  appropriate 
finite  difference  grids  are  illustrated  in  Figures  4  and  5.  Thus, the 
problem  of  evaluating  (3.12a  -  3.12d)  reduces  to  that  of  com¬ 
puting  the  xx  and  yy  components  of  the  electric  vector  potential 
dyads  integrated  over  the  aperture  pulses  ('  coordinates)  and 
matched  at  the  appropriate  pulse  centers.  The  major  obstacle  to 
these  computations  is  that  the  cavity  dyads  (Table  2)  converge 
more  and  more  slowly  as  the  distance  between  r_  and  r'  becomes 
small.  Summation  techniques  for  these  dyads  are  presented  in  the 
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next  section.  The  half-space  potential  terms  are  evaluated  by 
analytical  integration  of  the  singular  terms,  then  Gaussian 
quadrature  [Abramowitz  and  Stegun  1964,  P.  887]  is  used  to 
integrate  the  remaining  terms. 

Since  the  source  and  observation  points  in  (3.12e)  and 
(3.12f)  remain  distant  the  z  derivatives  may  be  brought  under 
the  integral  signs.  The  remaining  kernels  consist  of  exponen¬ 
tially  convergent  infinite  series  which  may  be  integrated  term 
by  term. 

The  derivatives  on  the  wire  potential  terms  (3.12g)  are 
kept  outside  the  integral  and  evaluated  by  finite  differences 
(Figure  6) .  The  major  obstacle  to  evaluation  of  this  term  is 
the  summation  of  a  cavity  potential  dyad  as  the  distance  between 
source  and  observation  points  becomes  small.  Techniques  for 
carrying  out  this  summation  are  presented  in  the  next  section. 


CHAPTER  4 


EVALUATION  OF  CAVITY  POTENTIAL  TERMS 


The  final  obstacle  to  the  solution  of  our  frequency 
domain  problem  is  the  matrix  fill  of  elements  involving  the 
cavity  vector  potentials  gF,  and  .  In  Table  2  these  dyads 
are  given  in  terms  of  their  threefold  infinite  eigenfunction 
expansions.  For  source  and  observation  points  distant  it  will 
be  seen  that  these  expansions  reduce  to  rapidly  exponentially 
convergent  double  series.  However,  as  source  and  observation 
points  become  close,  these  series  become  slowly  convergent. 

This  chapter  develops  an  efficient  means  of  evaluating  the 
potential  dyad  g^(r|r' )  for  aperture  source  and  observation 
points  close,  under  the  restriction  that  ka  (Figure  1)  is  less 
than  tt .  The  same  technique  will  be  used  to  evaluate  G^(r'  |r)  in 
(3.12g)  under  the  thin-wire  restriction. 

Consideration  is  now  given  to  the  xx  component  of  the 
dyad  gF  .  Similar  techniques  hold  for  g  921 .  From  Table  2 
with  z  »  z'  *  c  (Figure  1) 


(r  I  r ' )  = 


“  £}£ m£n 

£  —3 - k  sin (k  x) sin (k  x ' ) cos (k  y ) cos (k  v ' ) 

=n  -v*  xx  y  y 


KImn-* 

m=l 


(4.1a) 


S ( | x-x' | ,  y-y'S)  +  S(|x-x'|,  y+y')  -  S(x+x’,  y+y') 
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where 


S  (X,  Y) 


_ 1 

4at>c 


Z 

1=0 

m=l 

n=0 


ei,£m£n 


cos (k^X) 


cos(kyY)  . 


(4.1c) 


Since  the  terms  of  (4.1b)  depend  only  on  sums  and 
differences  of  the  respective  source  and  observation  points' 
components,  the  number  of  calculations  in  evaluating  (3.12b) 
and  (3.12d)  is  less  than  if  (4.1a)  were  applied  directly. 

For  source  and  observation  points  distant  the  threefold 
infinite  series  of  (4.1c)  may  be  reduced  to  a  double  sum  that 
yields  itself  to  rapid-  evaluation  by  digital  computer.  Equation 
(4.1c)  is  closed  in  Y  [Gradshteyn  and  Ryzhik  1965]  which  yields 


S  (X,  Y) 


®  e.  cosh  ( +  k2,  -  k2  (b-Y)] 

E  -■  ■  -  y .  —  ■  ■ — 

1=0  Al  +  k^  -  k2  sinh  +  k^  -  k2  b] 
x  y  x  y 

m=l 


cos (kxX) 
(4.2a) 


with  ka  (Figure  1)  less  than  it.  The  series  in  (4.2a)  is  exponen¬ 
tially  converging  and  readily  evaluated  provided  that  Y  is  not 
too  small  with  respect  to  the  aperture  dimensions  and  not  too 
close  to  2b.  Likewise  if  (4.1c)  is  closed  in  X  one  obtains 


S  (X,  Y)  - 


1 

4bc 


I 

1=0 

n=0 


ezen  T (X)  cos (k^Y) 


1 

4abc 


°®  e.e  cos(k  Y) 
-  I  n  y 

1=0  k“  +  k~  -  k2 
n=0  y  2 


(4.2b) 
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where 


cosh  [A2  +  k2  -  k2  ' (a-X)] 
_  y  2 

/kf  +  k^  -  k^  sinh  [  /k  +  k 
y  z  y  z 


Ta] 


if  k2  +  k2  -  k2  >  0 

y  z 


T  (X)  = 


-cos[/k2  -  k‘  -  k*  (a-X)] 


.2 

“JL 


T1 

‘z 


2  2 

kj  -  k‘  sin ( 


ky  -  *1  a» 


2  2 
if  k  ,  +  kz 

y  z 


2 

k  <  0 


Similar  techniques  may  be  used  to  close  the  second  term  of  (4.2b). 
The  result  is  a  single  finite  series  with  know  asympotic  form 
[Abramowitz  and  Stegun  1970,  p.  1005].  This  this  series  may 
readily  be  evaluated.  First  termwise  subtraction  of  the  corre¬ 
sponding  asymptotic  term  from  the  original  series  term  yields  a 
rapidly  converging  new  series.  Next  the  otherwise  slowly 
converging  asymptotic  series  is  added  to  the  modified  series  by 
its  known  analytical  form.  The  first  term  of  (4.2b)  is  exponen¬ 
tially  convergent  provided  X  is  not  too  small  or  too  close  to  2a. 

In  the  aperture  geometries  considered  (Figure  1) ,  the  only  case 
where  neither  (4.2a)  or  (4.2b)  yields  itself  readily  to  numerical 
evaluation  is  for  both  X  and  Y  small  compared  to  the  cavity 
dimensions.  Otherwise,  Seidel  [1977]  has  obtained  an  efficient 
method  for  numerical  evaluation  of  these  sums. 

The  remainder  of  this  chapter  is  concerned  with  evaluation 
of  (4.1c)  near  the  source,  that  for  X  *  |x-x' |  and  Y  =  |y-y' | 
small.  Equation  (4.1c)  is  closed  in  z,  making  use  of  Spiegel 
[1964,  p.  189]  and  Gradshteyn  and  Ryzhik  [1965,  p.  40].  This  yields 
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S{X/Y)  *  2ab  S  en  T(kx  +  ky  "  k2>  cos (k^X)  cosfk  Y)  (4.3) 

m=l  *  1 

n=o 

where 


This  series  in  (4.3)  has  an  asymptotic  form 


sASY<x,Y) 


.2  +  k2 

-X  +  ky 


cos(k  X)  cos(k  Y)  .  (4.4) 

x  y 


The  difference  series  whose  terms  consists  of  the  original  series 
terms  (4.3)  minus  the  corresponding  asymptotic  series  terms  is 
readily  evaluated  by  digital  computer.  Thus  (4.3)  may  be  evalu¬ 
ated  if  a  method  of  evaluating  the  otherwise  slowly  converging 
asymptotic  series  may  be  found. 


To  evaluate  SASY(X,Y) 


(4.4)  consider  the  single  series 


Sq(Y) 


(4.5) 


2  2  2.. 

where  Y  is  greater  than  zero  and  u  equals  kx  -  k  which  is 
greater  than  zero  under  the  restriction  that  ka  is  less  than  ir. 

An  alternate  expression  for  SQ(y)  will  be  obtained  from  Poisson's 


formula  [Papoulis  1962] 
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Equation  (4.8b)  follows  from  (4.8a)  plus  symmetry  considerations. 
The  integrals  of  (4.8b)  result  in  the  modified  Bessel  functions 
of  (4.8c)  [Morse  and  Feshbach  1953,  p.  1323].  Application  of 
(4.6a)  to  (4.7)  yields 


so(Y) 


-  I  £  [Ko(  |Y  +  2nb|A2  -  k2)  +  Ko(|2nb  -  Y|Vk2  -  k2)] 


(4.9a) 


=  ~  KQ(Vk2  ‘  k  Y)  +  IT  £  [K0([2nb  +  Y]A2  -  k2) 

n=l 


+  Kq( (2nb  -  Y] 


(4.9b) 


Applications  of  (4.9b)  to  (4.4)  yields 


sasy(x'y>  -  jtL  V^x  -  k5  Y>  cos(kxx) 

+  57  £  t  [K0i[2nb  +  Y]  -Al  -  k2) 

m=l  n=l 

+  K  ( [2nb  -  Y]  >4 2  -  k2) ]  cos (k  X).  (4.10) 

OX  X 

The  second  term  on  the  right  of  (4.10)  is  very  rapidly  convergent 
due  to  the  exponential  behavior  of  KQ's  asymptotic  expansion 
[Arfken  1971,  p.  517].  The  first  term  on  the  right  (4.10)  is 
exponentially  convergent  for  Y  >  0.  However,  as  Y  becomes  smaller 
the  series  converges  more  and  more  slowly.  A  summation  formula 
for  the  first  term  is  given  Lewin  [1951,  p.  83].  Thus 
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i?  E. 


k  Y) 


m=l 


cos (kxX) 


YQ(kY) 

4a 


+ 


1 

2tt 


cos [k/X2  +  Y2] 
/x2  +  Y2 


(4.11) 


costkr^2  +  (2na  +  x)2]  cos  [k/y2  +  (2na  -  X)2] 
A2  +  (2na  +  X)2  JF  +  (2na  -  X) 


Note  that  the  dominant  singularity  in  (4.11)  is  the  same  as  that 
of  the  half  space  dyadic  gF  of  Table  1.  Unfortunately  the 
infinite  series  (4.11)  is  slowly  convergent.  However  if  one 
approximates  each  term  of  this  series  by  a  finite  Taylor  series 
in  X  and  Y  about  the  origin,  the  summation  over  n  may  be  obtained 
by  analytical  methods.  For  example,  the  first  term  obtained  in 
the  Taylor  series  is  given  by 

V  C°-S^2-kan)  =  -  £n  [2  sin  (ka)  ]  (4.12) 

n=l  n 

[Abramowitz  and  Stegun  1970,  p.  1005].  Similar  formulas  may  be 
used  to  evaluate  the  higher  order  terms  in  the  Taylor  series 
approximation . 

Although  the  above  discussion  has  been  restricted  to  the 
evaluation  of  gF'  only  minor  changes  are  needed  to  evaluate 
gF,y^.  Also  the  terms  that  must  be  evaluated  (3.12)  are  integrals 
of  the  £&  and  ££  components  of  gF  .  The  sums  of  (4.2)  may  be 
evaluated  by  term  integration  [Titchmarsh  1939,  pp.  36-45].  The 
exponentially  convergent  terms  (4.10)  are  also  evaluated  termwise. 
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The  only  remaining  question  is  the  accuracy  of  the  previously 
discussed  Taylor  series  approximation  (4.11).  A  sixth  order 
Taylor  series  was  used  for  numerical  testing  and  yielded  no 
relative  error  worse  than  1*10  ^  in  these  terms  when  max(X,Y)/ 
a  _<  .05.  This  procedure  allows  rapid  evaluation  of  the  series 
on  the  right  of  (4.11),  provided  that  both  X  and  Y  are  small 
with  respect  to  the  cavity  dimensions  in  the  aperture  plane. 
Otherwise,  this  series  would  not  readily  yield  itself  to  numerical 
evaluation  on  a  digital  computer. 

The  above  summation  techniques  are  now  used  to  evaluate 
the  potential  terms  of  (3.12g).  Under  the  thin  wire 

approximation  all  current  is  assumed  to  reside  along  the  wire 
center.  Following  Lewin  [1975,  p.  121]  the  wire  radius  is  set 
to  zero  in  all  terms  except  where  needed  to  prevent  a  divergent 
series. 

To  reduce  the  number  of  computations  necessary  to 
evaluate  all  of  the  G^' ^  terms,  this  potential  is  expressed  as 
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This  series  is  rapidly  exponentially  convergent  provided  that 
Y  is  not  too  small  or  does  not  approach  2b.  Since 


S { Y)  =  S(2b-Y)  (4.16) 

all  that  remains  to  be  considered  is  y  small  with  respect  to 
the  cavity  dimensions. 

A  scheme  for  numerical  evaluation  of  S (Y)  with  y  small 

is  now  given.  Recall  that  F  =  (x  +r  ,y  ,  z  )  and  r'=  (x  ,y',z  ) 

n  c  w  ■‘m  c  c  c 

where  x  and  z  are  the  respective  x  and  z  coordinates  of  the 

C  G 

wire  center.  By  means  of  a  trig  identity  S(Y)  (4.14)  may  be 
expressed  as 


S(Y)  »  T2(Y)  -  Tl(Y) 


(4.17) 


where 


T1(Y) 


=  Y\  e.e  e.  - 

4abc  l  m  n  R2 


cos  [kx  (2xc+rw)  ]cos(k  Y)  (ss  )  z 


(4.18) 


m=l 

n=o 


Klmn  -  k 


and 


T2(Y) 


-Jr-  T 
4abG  i=l 
m=l 
n=o 


£„£  e.  cos  (k 
i  in  n  x 


rw)  cos(kyY)  (ss) 


(4.19 


K 


£mn 


-  k4 


Equation  (4.18)  may  be  closed  in  x  to  yield  a  rapid  exponentially 
converging  series.  The  wire  radius  in  (4.18)  is  set  to  zero. 

Next  (4.19)  is  closed  in  z  to  yield 


,  ®  e  e  cos (k  Y) sinh  fvk  +k  -k  (C-Z) . 

T  (V)  =  1  y  . -5L-S - 1 - * — 1 - 

2  4abc  i  7  9  /  9  9  9 

/k;  +  k"1  -  kz  sinh [ /k  +  k  -  kz  C]  (4.20) 

n=o  x  y  x  y 

This  series  has  the  asymptotic  from  SAgY/2.  (4.4).  Thus  the 
procedure  for  evaluation  of  this  term  has  already  been  given. 


CHAPTER  5 


SAMPLE  NUMERICAL  RESULTS 

The  numerical  solution  scheme  of  the  previous  sections 
was  implemented  by  means  of  a  fortran  program  on  a  Cyber  175 
computer.  In  all  examples  (Figures  7,10,14)  the  wire  radius  was 
set  at  one  millimeter,  the  number  of  aperture  pulses  in  the 
x-direction  was  4,  the  number  of  aperture  pulses  in  the 
y-direction  was  5,  and  the  number  of  wire  current  pulses  was  12; 
a  zero  half  pulse  counts  as  one  half. 

In  the  first  example  (Figure  7)  a  normally  incident 
plane  wave,  with  free- space  wavelength  equal  to  10  meters,  im¬ 
pinges  on  a  cavity  backed  aperture.  The  incident  electric 
field  is  y  directed  to  give  maximum  coupling  to  the  interior 
wire.  The  cavity  dimensions  are  a  =  1.4  meters,  b  =  4.5  meters, 
and  c  =  2.4  meters  in  the  x,  y,  and  z  dimensions  respectively. 

The  square  .02m  x  .02m  aperture  is  centered  with  midpoint 
r  =  (.7,2.25,1.2).  The  centered  thin  wire  runs  from  (.7,. 5, 1.2) 

to  (.7,  4., 1.2).  Figures  8  and  9  show  some  cuts  of  the  aperture 
electric  field.  As  expected  for  electrically  small  apertures 
the  aperture  electric  field  is  essentially  imaginary.  The 
real  part  is  seven  orders  of  magnitude  lower  than  the  imaginary 
part.  The  y  component  of  the  aperture  electric  field  has  even 
symmetry  about  both  axes  of  the  aperture,  while  the  x  component 
is  odd.  Figures  10  and  11  show  a  comparison  between  Seidel's 
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Figure  7  Problem  Geometry 


Figure  9.  Real  Aperture  Electric  Fields  For  Problem  Geometry 


Integral  Equation* 
Dipole  Moments 
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dipole  moment  solution  [Seidel  1977,  1978]  (dotted  line)  and 
the  integral  equation  approach  (solid  line).  Seidel’s  solution 
treated  the  aperture  coupling  by  the  use  of  dipole  moments  for 
elliptical  apertures,  in  this  case  circular  with  radius  .01  m  and 
located  inside  the  square  aperture.  Unfortunately  a  sign 
discrepancy  occurs  in  the  imaginary  part  of  the  current  (Figure 
11) .  Note  that  while  the  symmetry  of  the  integral 
equation  method  is  preserved  the  dipole  moment  has  lost  its 
symmetry.  At  first  one  might  suspect  that  the  imaginary  part 
of  the  current  (Figure  11)  is  in  the  noise,  however  after 
running  various  test  cases  this  sign  discrepancy  consistently 
occurred.  A  painstaking  search  of  our  computer  program  yielded 
no  error.  An  attempt  to  resolve  this  discrepancy  is  given  in 
Appendix  A. 

In  the  next  example  (Figure  12)  an  incident  plane 
wave  travels  in  the  y  direction,  is  polarized  with  a  z  directed 
incident  field,  and  has  free-space  wavelength  equal  to  one 
meter (lm) .  The  cavity  dimensions  are  .25m,  .3m,  and  .2m  in  the 
x,  y,  and  z  directions  respectively.  The  wire  is  oriented 
parallel  to  the  y  axis  and  runs  from  ( .  125 ,  .  05 ,  .  1)  to  (.125, 
.25,. 1).  The  square  aperture  (.02m  x  .02m)  is  centered  with 
midpoint  r  =  ( .125, .15,  .2) .  Figures  13  and  14  show  cuts  of 
the  aperture  electric  fields,  the  symmetry  of  the  aperture 
electric  fields  is  described  in  Table  3.  Figure  15  gives  a 
comparison  between  the  wire  currents  computed  by  Seidel's 


Figure  12  Problem  Geometry 


Figure  13.  Aperture  Electric  Fields  for  Problem  Geometry 
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Table  3  Aperture  Electric  Field's  Symmetry  for  Geometry  2. 


Component 

Real  or  Imaginary 

Aperture  x-axis 

X 

real 

odd 

even 

X 

imaginary 

odd 

odd 

y 

real 

even 

odd 

y 

imaginary 

even 

even 

Figure  15.  Wire  Currents  for  Problem  Geometry  2 
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dipole  moment  approach  and  the  technique  presented  in  this 

paper.  The  dipole  moment  solution  is  for  a  circular  aperture 

of  radius  .01m  and  with  centerpoint  r  =  ( .125, . 15 ,  .2)  . 

Again  the  sign  discrepancy  persists. 

Figure  16  has  the  same  geometry  as  Figure  7  except  the 

centered  aperture  is  made  rectangular  (.8m  x  ,3m).  A  study  of 

wire  current  versus  frequency  is  made  in  Figures  17  -  21. 

Figures  7,  8,  9,  and  20  show  the  variation  in  current  for 

4  4  8 

frequencies  of  10  .  4.10  -10  hertzs  respectively.  Note  the 

change  in  sign  of  the  wire  current  as  the  frequency  crosses 

7 

the  first  wire  resonance  of  4.16  •  10  hertzs.  This  resonant 
frequency  is  slightly  lower  than  that  for  a  wire  in  free  space 
Figures  21  and  22  show  the  variation  of  the  current  at  the 
wire  center  with  frequency.  Note  that  the  real  part  of  the 
current  (Figure  21)  is  proportional  to  the  frequency  squared 
in  the  quasi-static  portion  of  the  frequency  spectrum.  Also 
the  imaginary  part  of  the  current  (Figure  22)  is  proportional 
to  the  fifth  power  of  the  frequency  in  the  quasi-static  region 
The  amount  of  computation  time  for  each  frequency  point  is 
20  seconds  c.p.u.  time  on  the  Cyber  175  computer. 
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Figure  18.  Wire  Currents  for  Geometry  3. 


Figure  19.  Wire  Currents  for  Geometry  3. 


APPENDIX  A 


ON  THE  SIGN  DISCREPANCY 


In  this  appendix  analytical  techniques  are  combined 
with  numerical  results  in  an  attempt  to  resolve  the  sign 
discrepancy  of  Chapter  5.  In  particular  the  problem 
geometry  of  Figure  7  will  be  considered.  First  check  on 
the  signs  of  the  aperture  electric  fields  will  be  made  by 
comparison  with  currents  induced  when  a  normally  incident 
plane  wave  strikes  a  thin  perfectly  conducting  rectangular 
plate.  If  one  neglects  the  contributions  to  the  aperture 
fields  due  to  scattering  from  the  wire  and  cavity  walls 
these  problems  are  Babinet  equivalents. 

The  currents  induced  on  the  plate  (Figure  23) 


satisfy 
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where  r  is  located  on  the  plate  surface  P,  and  R  «  |r-r'|. 


(A. 2) 


For  a  normally  incident  polarized  plane  wave  (Figure  23) 
(A. 1)  and  (A. 2)  may  be  expressed  as 
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CA.4) 


If  one  neglects  the  contribution  of  the  aperture 
fileds  due  to  the  wire  and  side  and  back  cavity  walls  and 


specializes  the  incident  field  to  an  polarized  normally 


incident  field  (2.7a)  and  (2.7b)  become 
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1  P  v  C  P  P 

and  E^  =  H„  »  J.  where 


If  one  identifies  E 

y  2“y  2  - “X  _X  ~y 

2 

the  superscripts  "C"  and  "P"  denote  cavity  and  plate  respec- 


Txnc 


tively,  one  sees  that  for  identical  forcing  functions  =  1 

inc  C  IP 

and  E  =1,  the  aperture  E  should  be  -  y  • 

y  y  £  X 


Figure  24  shows  a  plate  scattering  problem  for  a  .15X 
inc 


square  plate  with  H  *  -1.  When  comparing  with  Figures  (8 


and  9)  with  E^nC  =  1  we  should  have  E^  to  have  the  same  sign 


as  Jx,  as  is  indeed  the  case.  Note  that  the  signs  of  the 


aperture  fields  and  wire  currents  should  not  change  in  the 
quasi-static  frequency  regime.  The  sign  discrepany  of  Chapter  5 
is  in  the  imaginary  part  of  the  wire  current.  Since  the 
Greens  functions  of  (2.7c)  are  real  the  imaginary  component  of 
the  wire  current  is  driven  by  the  real  part  of  the  aperture 
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field  Ey.  The  question  arises  as  to  whether  the  presence  of 
the  cavity  and  or  wire  changes  the  sign  of  the  real  part  of 
the  aperture  field  E^.  The  results  of  our  computations  indi¬ 
cate  that  this  is  not  the  case,  as  illusiurated  in  Figures 
(25-28) . 


The  imaginary  component  decreases  by  approximately 
factor  of  2  as  the  Green's  functions  kernel  in  CA.5)  and  (A. 6) 
changes  from 


1  ,cos(kR)  4  sin(kR)x  ,, 

_  <  -  l  g  )  IA.  n 

for  the  halfpace  to  approximately 

1  ,coskR  ..  sin(kR)  /7> 

2 F~  (  R  1  2R  n 

for  the  interior  cavity,  exterior  halfspace  problem.  It  is 
comforting  to  note  that  decreasing  the  imaginary  part  of  the 
kernel  also  reduced  the  real  part  of  E^,  since  an  entirely 
real  kernel  in  (A. 5)  and  (A. 6)  would  produce  a  zero  real 
component  of  both  Ey  and  Ex» 

An  argument  is  now  given  to  determine  the  sign  of  the 
imaginary  wire  current  (Figure  11).  Note  that  (2. 7C).  may  be 
expressed  as 


Jy(r)  =  iu)£  0  ( Ey )  (A.  7) 

where  0  is  a  real  linear  operator.  Note  that  the  real  and 
imaginary  parts  of  Ey  have  approximately  the  same  shape.  Also 
note  that  the  imaginary  part  of  Ey,  Ey  produces  an  negative 
real  current  (Figure  10).  Thus  0(Ey)  is  positive.  Since  the 


3 


63 


real  part  of  Ey(Ey)  has  approximately  the  same  shape  it  follows 
that  3(Ey)  is  again  positive,  so  that  the  imaginary  part  of 
Jy  (Figure  11)  should  be  positive. 
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